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INTERFEROMETRIC OPTICAL COMPONENT ANALYZER 
BASED ON ORTHOGONAL FILTERS 



BACKGROUND OF THE INVENTION 

Field of the Invention 

The principles of the present invention relate generally to the field of optical 
network analysis, and more particularly, to a system and method for determining the 
optical characteristics of an optical component under test. 

Description of Related Art 

An optical network analyzer is a vital tool for determining optical characteristics 
of optical components, such as fiber Bragg gratings. The optical characteristics 
determined by an optical network analyzer may include reflectance and transmittance of 
a particular two port or multiport optical component under test. The optical 
characteristics (e.g., reflection or transmission) of an optical component under test are 
described by a transfer function and are typically determined by measuring the 
amplitude and phase of optical signals that have been reflected by or transmitted 
through the component. The phase response characteristics of an optical component 
under test are often described by group delay or dispersion. 

Most conventional group delay measurement techniques were developed for 
optical network analyzers that utilize non-continuously tunable laser sources. A non- 
continuously tunable laser source is a laser source that can be tuned across a predefined 
range of frequencies in discrete frequency steps. A typical technique to measure 
amplitude, phase, and group delay is to apply an electrical stimulus in the form of 
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intensity or phase modulation to an optical signal, and then measure the electrical 
response of the transmitted or reflected optical signals using a well known phase 
sensitive electrical detection device, e.g., a lock-in amplifier or an electrical network 
analyzer. A concern with this technique is that the group delay measurement is indirect 
and limits the accuracy of the measurement. In addition, the technique requires a long 
measurement time to obtain an accurate result. Consequently, a long-term stability of 
the test setup is required to effectively utilize the technique. 

However, continuously tunable laser sources have recently become available. A 
continuously tunable laser source can continuously sweep a predefined range of 
frequencies without frequency jumps or mode hops. The availability of continuously 
tunable laser sources has allowed for development of interferometric methods for 
measuring the optical characteristics of optical components. The interferometric 
methods are based on direct measurements of phase differences between interfering 
optical signals. Typically, Fourier analysis of a heterodyne beat frequency directly 
related to the sweeping optical frequency of a continuously tunable laser source is used 
to measure the optical characteristics, including the group delay. A concern with the 
interferometric methods using Fourier analysis is that the frequency sweep of a 
continuously tunable laser source is non-uniform. The non-uniformity of the frequency 
sweep causes a similar non-uniformity in the resulting beat frequency that introduces an 
uncertainty in the calculation of the optical characteristics by means of the Fourier 
analysis. What is needed is a technique to determine the optical characteristics of the 
continuously tunable laser with minimal uncertainty. 
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SUMMARY OF THE INVENTION 

The principles of the present invention determine optical characteristics of an optical 
component. To overcome the problem of using tunable laser sources having non-uniform 
frequency sweeps to measure optical characteristics of optical components, an optical 
network analyzer or test system that utilizes amplitude and phase computing components, 
such as orthogonal filters, to compute the optical characteristics of an optical device under 
test (DUT) may be employed. The optical characteristics may include amplitude, phase, and 
group delay. The optical network analyzer may include an interferometer to measure the 
optical DUT, and a processor that processes the measured optical signals utilizing the 
orthogonal filters. The orthogonal filters include an in-phase and a quadrature filter. A 
reference interferometer may further be included to measure optical frequency and/or 
amplitude of the tunable laser. Amplitude and phase are computed from heterodyne beat 
signals produced by the interferometers, and group delay is computed based on the computed 
phase. 

One embodiment includes a system and method for measuring optical characteristics 
of an optical device under test (DUT). The system includes a light source for generating an 
optical signal applied to the optical DUT. A reference interferometer and a test 
interferometer are optically coupled to the light source. A computing unit is coupled to the 
interferometers, and utilizes the orthogonal filters in determining optical characteristics of the 
optical DUT. The optical characteristics may include a reflective transfer function, a 
transmissive transfer function, and group delay. By utilizing the principles of the present 
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invention, uncertainty in computing the optical characteristics of the optical DUT may be 
minimized. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram of an optical network analyzer configured to perform 
reflectance measurement in accordance with the principles of the present invention; 

Fig. 2 illustrates the frequency sweep of an ideal continuously tunable laser source 
with respect to time; 

Fig. 3 illustrates the AC coupled heterodyne beat signal when an ideal continuously 
tunable laser source is used; 

Fig. 4 illustrates a phase function of the heterodyne beat signal of Fig. 3; 

Fig. 5 illustrates the AC coupled heterodyne beat signal when a continuously tunable 
laser source that sweeps frequencies non-uniformly with respect to time is used; 

Fig. 6 illustrates a phase function of the heterodyne beat signal of Fig. 5; 

Fig. 7 is an exemplary block diagram for processing the optical signals measured by 
optical detectors according to Fig. 1 ; 

Fig. 8 illustrates a heterodyne beat signal as formed by the optical network analyzer of 

Fig. 1; 

Fig. 9 illustrates a filter function for in-phase detection of amplitude and phase of the 
heterodyne beat signal of Fig. 8; 
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Fig. 10 illustrates a filter function for quadrature detection of amplitude and phase of 
the heterodyne beat signal of Fig. 8; 

Fig. 1 1 is an exemplary block diagram for applying the filter functions to determine 
amplitude and phase of the heterodyne beat signal of Fig. 8 in the time domain; 
5 Figs. 12A-12D illustrate orthogonal filters in the time and frequency domains as 

5 applied in Fig. 11; 

yj Fig. 13 is another exemplary block diagram for computing amplitude and phase of the 

H= heterodyne beat signal of Fig. 8 in the frequency domain; 

^ Fig. 14 illustrates a transfer function of a positive frequency bandpass filter centered 

M at an angular frequency for computing amplitude and phase of the heterodyne beat signal of 
jtj Fi g- 8; 

sssi 

f=& Fig. 15 is another exemplary block diagram for computing amplitude and phase of the 

heterodyne beat signal of Fig. 8 in frequency domain using the transfer function of Fig. 14; 
Figs. 1 6 A and 1 6B illustrate graphs of in-phase and quadrature filters, respectively, 
15 used in computing amplitude and phase of the heterodyne beat signal of Fig. 8, where Fig. 
16B depicts a Hilbert transform; 

Fig. 17 is another exemplary block diagram for computing amplitude and phase of the 
heterodyne beat signal of Fig. 8 using the Hilbert transform of Fig 16B; 

Fig. 18 illustrates a positive frequency, all-pass filter transfer function utilized in 
20 computing amplitude and phase of the heterodyne beat signal of Fig. 8; 
Fig. 19 illustrates a time-to-frequency conversion function; 

6 

DALLAS2 777347v2 47429-00037 



PATENT APPLICATION 
PDN# 10010292-1 
47429-00037USPT 

Fig. 20 is an exemplary block diagram depicting a simple conversion from the time to 
frequency domain; 

Fig. 21 is an exemplary block diagram for computing a transfer function and group 
delay of the optical DUT utilizing a phase subtraction method; 
5 Fig. 22 is another exemplary block diagram for computing a transfer function and 

including group delay of the optical DUT utilizing a derivative ratio method; and 

|7j Fig. 23 is another exemplary block diagram for computing a transfer function and 

Co 

M= group delay of the optical DUT utilizing a clocking method. 

O 

Si 

il DETAILED DESCRIPTION OF THE EXEMPLARY EMBODIMENTS OF THE 
3 INVENTION 

fy 

Tunable laser sources as presently available frequency sweep in a non-uniform 
manner. A heterodyne beat signal having a beat frequency produced by an interferometer for 

15 measuring optical components is affected by the sweep non-uniformity of the tunable laser 
source. Typical variations of the heterodyne beat signal due to the sweep non-uniformity are 
substantially larger than those induced by the measured dispersion of the optical component. 
In characterizing optical components, typical optical measurement systems or optical network 
analyzers include two interferometers, a non-dispersive reference interferometer and a test 

20 interferometer that includes the optical component or device under test (DUT) being 

measured. The reference interferometer is used for measuring the non-uniform sweep of the 
tunable laser source, where the measurement of the non-uniform sweep is utilized for 
compensation of the measurement from the test interferometer. Heterodyne beat signals from 
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both interferometers may vary in frequency, phase, and amplitude. The precision of the 
amplitude and phase detection of the beat signals from the interferometers determines the 
accuracy of characterization of the optical component. 

In computing optical characteristics of the optical component, the heterodyne beat 
signals from the interferometers are utilized to compute optical frequency of the tunable laser 
source and a reflection or transmission transfer function, including amplitude and phase, of 
the optical DUT. In determining the amplitude and phase characteristics of the heterodyne 
beat signal, orthogonal filters may be utilized. The orthogonal filters are formed by in-phase 
and quadrature filters or other acceptable orthogonal filters in the time and/or frequency 
domain. Reflection or transmission transfer functions and group delay of the optical 
component may be determined from the amplitude and phase of the heterodyne beat signals. 
The reflection or transmission characteristics may be determined using various computational . 
methods, including subtraction, division, and clocking. 

With reference to Fig. 1, an optical network analyzer 102 in accordance with the 
present invention is shown. The analyzer measures optical characteristics of the optical 
device under test 104. In particular, the analyzer determines the amplitude of the reflection 
transfer function and the group delay of the optical DUT 104. As an example, the optical 
DUT may be a fiber Bragg grating, although the optical DUT may be any optical component. 
The optical network analyzer computes the amplitude and phase of the transfer function of 
the optical DUT in an efficient manner using an interferometric method based on orthogonal 
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filters of AC coupled heterodyne beat signals. The group delay is determined from the phase 
response. 

The optical network analyzer 102 includes a reference interferometer 106, a test 
interferometer 108, and a processing unit 110. As illustrated in Fig. 1, the reference and test 
interferometers 106 and 108 are configured as Michelson interferometers to measure the 
reflective characteristics of the optical DUT 104. However, other types of interferometers 
could instead be used. For example, the interferometers 106 and 108 may be configured as 
Mach-Zehnder interferometers to measure the transmissive characteristics of the optical DUT 
104. Alternatively, rather than using a reference interferometer, an optical or electronic 
counter as known in the art may be utilized to detect the phase of the continuously tunable 
laser source. The reference and the test interferometers 1 06 and 1 08 both utilize a 
continuously tunable laser source 112. The continuously tunable laser source 1 12 is 
configured to generate an optical signal having a frequency that continuously sweeps across a 
predefined range of frequencies in a non-uniform manner. That is, the rate of change in 
frequency of the optical signal generated by the tunable laser source varies with time as the 
predefined range of frequencies is swept. 

As explained in detail below, the transfer function of the optical DUT 104 is 
computed by measuring the amplitude and the phase changes of an AC coupled heterodyne 
beat signal from the test interferometer 108 caused by the optical DUT 104. The heterodyne 
beat signal is produced from an input optical signal generated by the continuously tunable 
laser source 112. However, due to the frequency sweep non-uniformity of the laser source, 
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the optical frequency changes of the input optical signal are non-uniform. The non- 
uniformity introduces undesired frequency/phase changes into the heterodyne beat signal. 
The reference interferometer 106 and the processing unit 110 form an optical frequency 
counter that measures the undesired phase changes caused by the non-uniform optical 
frequency sweep of the input optical signal. Alternative optical frequency counters as known 
in the art may be utilized. The reference interferometer 106 produces a reference heterodyne 
beat signal that includes the phase changes caused by the frequency sweep non-uniformity of 
the tunable laser source. Additionally, the reference heterodyne beat signal may include 
amplitude caused by the tunable laser source. The reference heterodyne beat signal is then 
detected to measure the undesired phase changes. Since the measured phase changes are 

mi 

jy common to both interferometers 106 and 108, the distortions due to the non-uniform 
jU» frequency sweep can be removed from the heterodyne beat signal produced by the test 

interferometer 108 using the reference heterodyne beat signal from the reference 

interferometer 106. 

15 The reference interferometer 106 of the optical network analyzer 102 includes the 

continuously tunable laser source 1 12, an optical coupler 114, Faraday mirrors 116 and 118, 
and a reference detector 120. These components of the reference interferometer 106 are 
connected to each other by optical fibers 122, 124, 126, and 128. The optical fiber 122 
connects the tunable laser source to the optical coupler 1 14, while the optical fiber 126 

20 connects the reference detector 120 to the optical coupler 1 14. Similarly, the optical fiber 
124 connects the optical coupler 1 14 to the Faraday mirror 118, while the optical fiber 128 
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connects the optical coupler 1 14 to the Faraday mirror 116. The optical coupler 1 14 and the 
optical fibers 122, 124, 126, and 128 are single mode (SM) components. The length of the 
optical fiber 128 differs from the length of the optical fiber 124 by the distance AL* that 
introduces delay into the reference interferometer 106 and defines its free spectral range 
(FSR). 

In operation, the tunable laser source 112 continuously generates an optical signal in a 
sweeping range of frequencies. As an example, the tunable laser source 112 may frequency 
sweep over an approximately 140 nm wavelength range starting from a wavelength of 1470 
nm or 1510 nm wavelength. The optical signal generated by the tunable laser source 1 12 is 
transmitted to the optical coupler 114 through the optical fiber 122. The optical coupler 114 
then divides the optical signal into two optical signals so that the divided optical signals are 
transmitted to the Faraday mirrors 116 and 118 through the optical fibers 128 and 124, 
respectively. The optical signal transmitted through the optical fiber 128 is reflected back to 
the optical coupler 1 14 by the Faraday mirror 1 16 at the end of the optical fiber 128. 
Similarly, the optical signal transmitted through the optical fiber 124 is reflected back to the 
optical coupler 1 14 by the Faraday mirror 1 18 at the end of the optical fiber 124. The 
Faraday mirrors 116 and 118 reflect light in respective substantially orthogonal polarization 
states, thereby ensuring that the returning optical signals have nearly the same polarization 
state. 

Due to the extra distance traveled by the optical signal in the optical fiber 128, the 
frequency of the returning optical signal in the optical fiber 128 differs from the frequency of 
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the returning optical signal in the optical fiber 124 at a given moment in time. As illustrated 
in Fig. 1, the frequency of the returning optical signal in the optical fiber 124 is denoted by 
vi, and the frequency of the returning optical signal in the optical fiber 128 is denoted by v 2 . 
Assuming the tunable laser source 112 sweeps the predefined range of frequencies at the rate 
Y, the optical frequencies vi and v 2 may be expressed as: 

v/ = yt + v 0 

V2 = yt - yx + vo , 

where v G is the initial optical frequency of the tunable laser source, and x is the round trip 
delay in the length AL R of the optical fiber 128. The round trip delay x is defined as: 
2nAL R 

T , 

C 

where n is the refractive index of the optical fiber 128 and c is the speed of light in a vacuum. 
Assuming the constant sweep rate, y, the frequency difference between the returning optical 
signals in the optical fibers 124 and 128 is: 

Av = yr 

The reference interferometer 106 is assumed to be non-dispersive, i.e., r does not 
depend on the optical frequency of the optical signals. Alternatively, the reference 
interferometer 106 may be dispersion compensated. The phase difference between the two 
reflected signals may then be expressed as: 

0 = 2/rA vt + 0 O = iTtyrt + 0 O , 
where 2tcAv is the angular frequency difference. 
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The returning optical signals in the optical fibers 124 and 128 are combined at the 
optical coupler 114 and transmitted to the reference detector 120 through the optical fiber 
126. When combined, the optical signal returning from the optical fiber 124 interferes with 
the optical signal returning from the optical fiber 128. The intensity of the interfering optical 
signals observed at the reference detector may be expressed as: 

/ = I a + I Q cos(2/ryrt + 0 O ) 
with the assumption that the polarization state of the returning optical signals is the same. 
The term I 0 cos(2^ryxt + 0 O ) defines the heterodyne beat signal produced by the returning 
optical signals. At the reference detector 120, the heterodyne beat signal is AC coupled by 
removing the DC component of the signal. The AC coupled heterodyne beat signal is 
transmitted to the processing unit 1 10 for signal processing (i.e., determining group delay of 
the optical DUT 104). 

The test interferometer 108 of the optical network analyzer 102 includes the 
continuously tunable laser source 1 12, an optical coupler 130, a mirror 132, the optical DUT 
104, and a DUT detector 134. Similar to the reference interferometer 106, the components of 
the test interferometer 108 are connected to each other by optical fibers 136, 138, 140 and 
142. The optical fiber 140 connects the tunable laser source 1 12 to the optical coupler 130, 
while the optical fiber 136 connects the DUT detector 134 to the optical coupler 130. 
Similarly, the optical fiber 142 connects the optical coupler 130 to the mirror 132 that may be 
a gold mirror, while the optical fiber 138 connects the optical coupler to the optical DUT 104. 
An additional polarization controller may be used in one of the arms to match the polarization 
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states of the returning optical signals. The optical coupler 130 and the optical fibers 136, 138, 
140 and 142 are also single mode (SM) components. The length of the optical fiber 138 
differs from the length of the optical fiber 142 by the distance AL& that introduces delay into 
the test interferometer 108 and defines its free spectral range. 
5 The operation of the test interferometer 108 is similar to the operation of the reference 

*£l interferometer 106. However, the round trip delay xt for the test interferometer 108 is 
iy dependent on the frequency of the optical signal from the tunable laser source 112 due to the 
dispersion within the optical DUT 104, or expressed mathematically, x T = tt(v). In addition, 
^ the amplitude of the heterodyne beat signal depends on transmission or reflection properties 

U of the DUT. The amplitude and phase response of the optical DUT 104 may be derived by 

fy 

fy measuring changes in phase and amplitude of the heterodyne beat signal observed at the DUT 
jM= detector 134. 

The phase of the heterodyne beat signal of the test interferometer 108 is modified by 
the sweep non-uniformity of the tunable laser source 112. That is, the fact that the tunable 

1 5 laser source sweeps frequencies in a non-uniform manner with respect to time introduces 
additional phase shift into the phase of the heterodyne beat signal detected at the DUT 
detector 134 of the test interferometer 108. The same phenomenon exists for the reference 
interferometer 106. Therefore, both interferometers contain information about the non- 
uniform frequency sweep of the tunable laser source 112. The effect of the non-uniform 

20 frequency sweep by the tunable laser source is described below with reference to Figs, 2-6, 
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Assuming that the tunable laser source 1 12 is an ideal tunable laser source that 
continuously sweeps the predefined range of frequencies (e.g., optical C- and L-bands) 
uniformly with respect to time, the change of frequency with respect to time for the optical 
signals generated by the laser source is linear, as illustrated in Fig. 2. Consequently, if the 
5 device under test has no dispersion, then the resulting AC coupled heterodyne beat signal is a 
5 pure sine wave of constant frequency as shown in Fig. 3. The phase, cp(t), of the heterodyne 

J53. 

IM beat signal is linear as illustrated in Fig. 4. 

=Z However, for the tunable laser source 112 that sweeps frequencies non-uniformly, the 

w frequency of the frequency of the heterodyne beat frequency varies, as illustrated in Fig. 5. 
i§ Consequently, the phase cp(0 is non-linear, as illustrated in Fig. 6. Therefore, changes in 
fy phase of the heterodyne beat signal at the DUT detector 134 and at the reference detector 120 
H are strongly influenced by the sweep non-uniformity of the tunable laser source. The phase 
changes due to the dispersive characteristics of the optical DUT 104 observed at the DUT 
detector 134 are indistinguishable from and dominated by phase changes induced by the 
1 5 sweep non-uniformity of the tunable laser source 112. 

The optical network analyzer 102 resolves the problem of phase changes being 
generated by the non-uniform sweep frequency of the tunable laser source 1 12 by 
determining phase of the heterodyne beat signal from the reference interferometer 106 to 
effectively "compensate" for the non-uniformity of the sweep of the tunable laser source 1 12. 
20 This compensation leads to the recovery of the true linear optical frequency scale as 
performed in classical network analyzers. The quality of the optical network analyzer 
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depends on the precision of the recovered amplitude and phase of the heterodyne beat signal 
s(t). 

Four exemplary orthogonal filter embodiments are provided for computing amplitude 
and phase of the heterodyne beat signals s(t) from the reference and test interferometers 106 
5 and 108. The four embodiments of the orthogonal filters include: (i) in-phase and quadrature 
sj3 filtering in the time domain utilizing convolution, (ii) in-phase and quadrature filtering in the 
frequency domain, (iii) single sided filtering, and (iv) all-pass filtering utilizing Hilbert 
transforms. 

^ Fig. 7 is an exemplary block diagram 700 for measuring and processing the 

P$i heterodyne beat signals s(t) produced by the reference and test interferometers 106 and 108 of 

|*hJ 

fG Fig- 1- The optical fibers 126 and 136 are optically coupled to reference 120 and DUT 134 

k& detectors, respectively. The reference 120 and DUT 134 detectors are coupled to 

preamplifiers 705a and 705b. The preamplifiers are further coupled to analog-to-digital 
(A/D) converters 710a and 710b that are also coupled to the processing unit 110. The 

15 processing unit 110 may include a reference processor 715a to process data from the 

reference detector 120 and a test processor 715b to process data from the DUT detector 134. 
Each of the processors 715a and 715b are coupled to a common processor 715c. The 
processors 715a, 715b, and 715c maybe coupled to internal memory 720a, 720b, and 720c 
(cumulatively 720), respectively. Alternatively, the memory 720 may be external from the 

20 processors 715 a- 7 1 5c. 
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In operation, the reference and DUT detectors 120 and 134 receive optical signals 
from the reference and test interferometers 106 and 108 (not shown) and convert the optical 
signals into electrical signals. The electrical signals may represent heterodyne beat signals 
generated by the interferometers 106 and 108. The electrical signals are communicated to the 
processing unit 1 10 for processing the heterodyne beat signals and computing the transfer 
function and group delay of the optical DUT 104. 

The preamplifiers 705a and 705b receive the electrical signals from the reference and 
DUT detectors 120 and 134, respectively. The preamplifiers 705a and 705b may be used to 
filter and/or scale the electrical signals. The A/D converters 710a and 710b convert the 
analog signals received from the preamplifiers 705a and 705b into digital signals. The digital 
signals are received from the A/D converters 710a and 710b by the processors 715a and 
715b, respectively, to compute amplitude and phase of the heterodyne beat signals produced 
by the reference and test interferometers 106 and 108. The amplitude and phase signals may 
be computed by utilizing orthogonal filters, and are received by the common processor 715c 
that computes the transfer function and group delay of the optical DUT 104. It should be 
understood that the processing unit 110 may alternatively be a single processor for 
performing the computations. The processors 715a-715c may be general purpose processors 
or specialized processors, such as digital signal processors (DSP). Alternatively, digital 
devices other than processors, such as a programmable logic gate array, may be utilized in the 
processing unit 1 10 to perform the computations. In another embodiment, rather than 
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converting the analog signals to digital signals, the processing unit 110 may determine the 
transfer function and using analog components. 

To further understand the heterodyne beat signals s(t) received from the 
interferometers 106 and 108 and how the heterodyne beat signals are utilized in processing 
the transfer function and of the optical DUT 104, the following description is provided. A 
frequency difference between the interfering optical waves giving rise to the heterodyne beat 
signal is Av = / B = yT T , where y is the optical frequency sweep rate and x T is the round trip 
delay of the test interferometer 108. For the reference interferometer 106, t r is independent 
of the optical frequency (i.e., t r (v) = x R = constant). For the test interferometer 108, x T may 
change with the optical frequency produced by the tunable laser source 1 12 to provide 
information about dispersive properties of the optical DUT 104. In general, the sweep rate y 
is not constant (i.e., y = y(t)). Therefore, the heterodyne beat frequency changes with time 
and optical frequency. The heterodyne beat frequency may be computed as: 

t-r T {v) 

The reference and test interferometers 106 and 108 may be used to distinguish 
between the phase changes in the heterodyne beat frequency due to the non-uniform sweep of 
the tunable laser source 1 12 and dispersion of the optical DUT 104. The heterodyne beat 
signal measured by the test interferometer 108 varies not only in frequency, but also in 
amplitude as reflectance or transmittance of an optical DUT 104 changes with the optical 
frequency. Therefore, the signals measured have an oscillatory character with a modulated 
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amplitude and frequency (phase). The precision with which the amplitude and the phase of 
the heterodyne beat signals is detected establishes performance of the optical component 
analyzer. 

Fig. 8 illustrates an exemplary heterodyne beat signal s(t) 800 from the DUT 
5 interferometer with a varying amplitude and frequency. The varying amplitude, shown by 
S envelope 805, has a minimum amplitude ai and a maximum amplitude a 2 , where ai < a 2 . The 
W varying frequency of the heterodyne beat signal s(t) has a maximum frequency fi and a 

minimum frequency f 2 , where fi > f 2 . 
O Orthogonal filters may be used to measure the heterodyne beat signal s(t) of Fig. 8. In 

§ particular, orthogonal filters can be used to recover the amplitude and phase of the 
111 heterodyne beat signal s(t). A spectral component of a signal at a frequency fb can be 



determined by calculating a correlation of the signal at frequency f 0 with a signal periodically 
varying at the frequency f 0 (i.e., with cos(27if 0 t)). The spectral component may be computed 
as: 



The broader spectrum of frequencies may be detected by selecting a periodically varying 
function that does not extend from -oo to +oo, but rather is limited by an envelope, h(t), 
decaying to zero for final values of time. 

The first exemplary embodiment of the orthogonal filters includes in-phase and 
20 quadrature filters in the time domain as provided in FIGS. 9 and 10. Fig. 9 illustrates an in- 



15 




DALLAS2 777347v2 47429-00037 



PATENT APPLICATION 
PDN# 10010292-1 
47429-00037USPT 



phase filter for detecting a frequency band center around fo. The function in Fig. 9 is 
described by p(t) = h(t)cos(27tf 0 t). The envelope 905 is described by a function h(t). Because 
correlation of the heterodyne beat signal with the function p(t) recovers the original 
heterodyne beat signal s(t) without changing its phase, the filter function p(t) describes the in- 
5 phase filter. The in-phase component of the heterodyne beat signal may therefore be 
.yj computed by: 

jjy x(t)= \s(t + T)p(T)dr . 



Fig. 10 illustrates a quadrature filter that may be used for the detection of a frequency 
f3 band centered around fb. The quadrature filter is defined in a similar manner as the in-phase 
m filter. The phase shift of n/2 is realized by replacing the cosine function of the in-phase filter 
C with a sine function. The envelope h(t) is the same as the envelope 905 of the in-phase filter, 
thereby maintaining the same spectral response of the in-phase filter 900. The equation q(t) = 
h(t)sin(27tf 0 t) describes the quadrature filter. Correlation of the heterodyne beat signal s(t) 
with the function q(t) shifts the phase of the heterodyne beat signal s(t) by nil: 

15 y(t) = ]s(t + T)q(T)dr , 

-co 

where y(t) is the quadrature component of the heterodyne beat signal s(t). 

Once the in-phase and quadrature components are known, then the amplitude and phase of 

the incoming signal can be calculated as: 



r(t) = Jjc 2 (0 + /(0 (amplitude) 

20 and 
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<p{t) = arctan^^^J (phase). 



It is worth noting that for the selected in-phase p(t) and quadrature q(t) filters, 
correlation is the same as convolution with exception of the sign of q(t). Since sign is 
irrelevant in the instant process of determining the optical characteristics, either correlation or 
.£§ convolution may be utilized. 

nil Fig. 1 1 is an exemplary block diagram for computing amplitude and phase of the 

;: m Z 

S3 heterodyne beat signal s(t) in time domain. The in-phase filter p(t) is convolved with the 
2 heterodyne beat signal s(t) at 1 105a to produce a time domain signal x(t), where the time 
:L% domain signal is the in-phase component of the heterodyne beat signal. The quadrature filter 
IS q(t) is convolved with the heterodyne beat signal s(t) at 1 105b to produce a time domain 
5 signal y(t), where y(t) is the quadrature component of the heterodyne beat signal s(t). The 

quadrature signals, x(t) and y(t), are utilized by function 1 1 10 to produce the amplitude and 

phase of the heterodyne beat signal s(t). 

In the time domain, the in-phase filter p(t) is mathematically represented as a product 
1 5 of the envelope function h(t) and cos(27tf 0 t). 

p(t) = h(t) cos(2itfct) = h(t) cos(G> 0 t), 

where coo is the angular frequency. The above equation can be rewritten as: 



p{t) = hit) 



e -Jo> 0 t + e j<o 0 t 



21 
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The second exemplary embodiment of the orthogonal filters includes in-phase and 
quadrature filters in the frequency domain. Assume that the Fourier transform of h(t) is 
represented by H(co), i.e., F[h(t)] = H(oo), then 



The function H(co - a>o) represents a bandpass filter centered at the angular frequency 
cd 0 . The shape of the band-pass filter is uniquely defined by the envelope function h(t) in the 
time domain representation of the filter. It should be understood that the narrower the 
envelope h(t) is in the time domain, the wider the filter H(co) is in the frequency domain. 
Because of this time- frequency relationship, the selection of the quickly decaying envelope 
h(t) in Figs. 9-10 results in the wide passband in the frequency domain. 

An analogous procedure is used to derive the quadrature filter Q(co), wherein: 



The quadrature filter Q(co) provides the same frequency response H((o-(o 0 ), however, the 
phase is shifted by nj2 . Since in-phase filter p(t) is real and even, the Fourier transform 
P(co) is also real and even. The quadrature filter q(t) is real and odd, therefore, the Fourier 
transform Q(co) is odd and imaginary. 

Figures 12A and 12B are the time and frequency domain representations of the in- 

— i r: u^+u +v»^ +;*v% A <**A fi-omiPtimrHnmciinc r^cn^r.ti vp.1v Ficrft '\').C*. and 1 7 T) Henint 

pllCl&C 111 U^J til CUlVt llV^JUVilVJl v*vAii***** — 5 • O-* -*- *- 1 

the quadrature filter in the time and frequency domains, respectively. 



P(») = F[p(t)] = M±KQe- J '» +^h(t)e 
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When implementing the orthogonal filters in the frequency domain, the heterodyne 
beat signal s(t) is: (i) transformed to the frequency domain, (ii) filtered using P(co) and Q(co), 
and (iii) transformed back to the time domain. One reason for implementing the orthogonal 
filters in the frequency domain is that the computations in the frequency domain, which 
include addition and multiplication, is computationally faster and more efficient than 
computations in the time domain, which include convolution. The time-to-frequency and 
frequency-to-time conversions are performed because convolution may be performed as 

tl multiplication in the frequency domain. This process is illustrated in FIG 13 in form of a 

O 

" block diagram 1 300 and constitutes another exemplary embodiment for computing the 
H amplitude and phase of the heterodyne beat signal s(t). The heterodyne beat signal s(t) is 

transformed by the Fourier transform 1305 into the frequency domain S(co). The heterodyne 
^ beat signal S(co) in the frequency domain is filtered by an in-phase filter P(co) 1310a and a 
quadrature phase filter Q(co) 1310b in order to generate signals X(co) and Y(co), respectively. 
The signals X(a>) and Y(a>) are transformed from the frequency domain back to the time 
1 5 domain via an inverse Fourier transform 1315a and 1315b. The time domain signals x(t) and 
y(t) resulting from the inverse Fourier transform operations are received by the function 
1110. The function 1110 computes the time domain amplitude r(t) and phase cp(t) of the 
heterodyne beat signal s(t). 

The third exemplary embodiment of the orthogonal filters includes single-sided 
20 filtering. The orthogonal filters may be represented by a transfer function T(co) having real 
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and imaginary parts defined by P(co) and Q(co), respectively. Mathematically, the transfer 
function T(a>) may be described as 
T(©) = P(©)-jQ(fl>). 

Output z(t) of the single-sided filter, in response to the heterodyne beat signal s(t), is 
J> imaginary. The signal x(t) is the real part of the output z(t), while the signal y(t) is the 
^ imaginary part of z(t), where z(t) may be expressed as z(t) = x(t) - jy(t). The signals x(t) and 
m y(t) 5 are in quadrature (i.e., 90 degrees out of phase). 

q The output z(t) is an analytic signal with respect to x(t), and may be expressed as 

-L. z(t) = r(t)exp(jcp(t)X 

CO , yft\ 

M where the amplitude r(t) =^x 2 (t) + y 2 (t) and the phase <p(t) = arctan^- . Typically, the 

sign of the phase cp(t) is irrelevant. The amplitude and phase of the heterodyne beat signal' 
s(t) may be found as the absolute value of z(t) and the argument of z(t), expressed as 
r(t) = |z(0| and <p(0 = arg(z(f)) . By using expressions for P((o) and Q(co), the transfer 
function T(ca) may be expressed in terms of function H(<o). Therefore, the transfer function 

1 5 T(co) is uniquely described by H(co-co G ) and, as before, represented as a band-pass filter being 
centered at the angular frequency co 0 . However, the transfer function or filter T(co) operates 
on the positive frequencies only while the negative frequencies are suppressed (i.e., single- 
sided filtering). In the time domain, the function t(t) is complex with its real and imaginary 
parts being described by the in-phase p(t) and quadrature q(t) filters. The transfer function 

20 T(co) is illustrated in Fig. 14. 
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Fig. 15 is another exemplary block diagram for computing the amplitude and phase of 
the heterodyne beat signal s(t). The heterodyne beat signal s(t) is transformed into the 
frequency domain, S(co) by the Fourier transform 1305. The heterodyne beat signal S(<o) is 
received and filtered by a transfer function T(co) 1505 to produce signal Z(cd). The signal 
5 Z(co) is transformed to the time domain via the inverse Fourier transform 1315. The function 
3 1110 produces the amplitude r(t) and phase (p(t) of the heterodyne beat signal s(t). 
|7j The fourth exemplary embodiment of the orthogonal filters includes performing all- 

.J3J-A 

H pass filtering. The all-pass filters are a special case of orthogonal filters. The sweep rate of 
e the tunable laser typically varies around some nominal sweep rate y, resulting in the 
W heterodyne beat frequency varying around yx T , where xt is a round trip delay of the test 
m interferometer 108. Therefore, the natural selection of the filter type is a band-pass filter. 
Ju However, an all-pass filter may also be considered. In implementing an all-pass filter, the 

quadrature filter q(t) may be implemented by a Hilbert transform. Shown in Figures 16A and 
16B are all-pass orthogonal filters in the frequency domain, where the in-phase and 
15 quadrature filters are illustrated. As shown, the in-phase filter P(co) is real and even. The 
quadrature filter Q(a>) is shown to be imaginary and odd. The exemplary block diagram for 
the all-pass orthogonal filters in the frequency domain is the same as that in FIG. 13. As 
before, the heterodyne beat signal s(t) remains unchanged after being processed by the in- 
phase filter. The phase of the heterodyne beat signal s(t) is shifted by n/2 for all frequencies 
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by the quadrature filter Q(co). The resulting signal y(t) may be found in the time domain by 



Fig. 17 is another exemplary block diagram for computing the amplitude and phase of 
the heterodyne beat signal s(t) utilizing the Hilbert transform. As shown, the heterodyne beat 
signal is received directly by the function 1110 and via the Hilbert transform 1405. The 
Hilbert transform 1405 performs the quadrature filtering on the heterodyne beat signal s(t) so 
that the function 1110 may compute the amplitude and phase of the heterodyne beat signal 



The all-pass filter transfer function T(co) takes a simple form. The all-pass filter 
suppresses the negative frequencies, while leaving the positive frequencies unchanged. This 
can be seen in Fig. 18 by the transfer function being shown as having zero for the negative 
frequencies and unity for the positive frequencies. 

To this point, several methods for computing (i.e., recovering) the amplitude and 
phase of the heterodyne beat signal s(t) using the four exemplary embodiments of orthogonal 
filters have been shown, including (i) in-phase and quadrature filtering in the time domain 
utilizing convolution, (ii) in-phase and quadrature filtering in the frequency domain, (iii) 
single sided filtering, and (iv) all-pass filtering utilizing Hilbert transforms. Once the 
amplitude and phase information of the heterodyne beat signals are computed by utilizing 
orthogonal filters (of any form) from both interferometers 106 and 108, the optical 
characteristics of the device under test may then be computed. As shown in Fig. 1, both the 




s(t). 



26 



DALLAS2 777347v2 47429-00037 



PATENT APPLICATION 
PDN# 10010292-1 
47429-00037USPT 



reference and test interferometers 106 and 108 are used to generate amplitude and phase 
information. The generated amplitude and phase signals are: r R (t) (amplitude of the 
heterodyne beat signal from the reference interferometer 106), (pR(t) (phase of the heterodyne 
beat signal from the reference interferometer 106), r T (t) (amplitude of the heterodyne beat 
5 signal from the test interferometer 108), and <p T (t) (phase of the heterodyne beat signal from 
jyS the test interferometer 108). 

ty Since the reference interferometer 106 is non-dispersive or compensated for 

^ dispersion, the recovered phase is proportional to the optical frequency. Every 2n radians of 

S3 

' u phase change corresponds to the change in the optical frequency equal to the free spectral 

H range (FSR) of the interferometer. Therefore, the function 9 R (t) provides a unique 

|y conversion between time and the optical frequency. 

where v is the optical frequency, and x R is the delay in the reference interferometer 106. 

Fig. 19 is an exemplary conversion curve or function 1900 for converting between 
15 time and the optical frequency. The conversion curve 1900 is non-uniform and describes the 
non-uniformity of the tunable laser source. Fig. 20 depicts a simple block diagram of a 
process 2000 for a conversion from the time to the optical frequency domain. 

Conversion from the time to the optical frequency domain may be performed by a 
wavelength counting device. In particular, the conversion may be performed by the reference 
20 interferometer. The heterodyne beat signal s(t) from the reference interferometer 106 may 
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also be used for clocking the acquisition of the heterodyne beat signal from the test 
interferometer 108. In this embodiment, the heterodyne beat signal from the test 
interferometer is sampled in equal optical frequency increments. Hence, the signals 
recovered by the orthogonal filters depend on the optical frequency v. The recovered signals 
5 r(v), and (p(v) are the amplitude and phase, respectively, of the heterodyne beat signal from 
J3 the test interferometer 106. 

y Three exemplary methods for determining optical characteristics of the device under 

test, including: (i) amplitude r(v), (ii) phase 0(v), and (iii) group delay x g (v), are shown that 
T utilize the heterodyne beat signals s(t) produced by the reference and test interferometers 106 
CO and 108. These methods include (i) phase subtraction, (ii) derivative ratio (i.e., division), and 
CU (iii) clocking. 

N 1 Fig. 21 shows an exemplary block diagram 2100 that utilizes the phase subtraction 

method for computing the group delay of the optical DUT 104. As shown, the phase cp R (t) 
generated by the reference interferometer 106 is weighted by a=T T /T R and subtracted from the 

15 test phase cp T generated by the test interferometer 108. The subtraction removes any sweep 
non-uniformity related to phase perturbations of the tunable laser source. Consequently, the 
phase 0(t) substantially contains information about dispersion. The phase 0(t) is converted 
to an optical frequency function 0(v) and is then differentiated to find the group delay t g (v). 
The amplitude r T (t) is also converted to a function r T (v) in the optical frequency domain. 

20 Fig. 22 is an exemplary block diagram 2200 that utilizes the derivative ratio method 

for computing the group delay of the optical DUT 104. The phase functions cp T (t) and (p R (t) 
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are differentiated by differentiators 2205a and 2205b, respectively, and divided by divider 
2210. The process of division removes sweep non-uniformity related perturbations. The 
resulting function is equal to the group delay x g (t) after multiplication x R . The group delay 
Tg(t) and amplitude r T (t) are converted from the time to the optical frequency functions by the 
conversion function 2000. 

Fig. 23 is an exemplary block diagram 2300 that utilizes the clocking method for 
computing the group delay x g (v) of the optical DUT 104. The heterodyne beat signal s(t) of 
the reference interferometer may be used to clock the acquisition of the heterodyne beat 
signal from the test interferometer 108. Such an implementation leads to the signal being 
sampled uniformly in the optical frequency domain, v. The orthogonal filters are used to 
recover the optical frequency functions <p T (v) and r T (v). The group delay x g (v) is found by 

computing the derivative 2015, d( Pr {y ^ 9 of the phase of the heterodyne beat signal from the 

dv 

test interferometer 108. 

The previous description is of exemplary embodiments for implementing the 
principles of the present invention, and the scope of the invention should not necessarily be 
limited by this description. The scope of the present invention is instead defined by the 
following claims. 
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